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ABSTRACT 

We define a polynomial measure of multiparticle entanglement which is scalable, i.e., which 
applies to any number of spin-| particles. By evaluating it for three particle states, for 
eigenstates of the one dimensional Heisenberg antiferromagnet and on quantum error cor- 
recting code subspaces, we illustrate the extent to which it quantifies global entanglement. 
We also apply it to track the evolution of entanglement during a quantum computation. 
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Although entanglement has been recognized as a remarkable feature of quantum me- 
chanics since Schrodinger introduced the word [1] in response to Einstein, Podolsky and 
Rosen's famous paper [2], it remains only incompletely understood. In fact, for more than 
two particles — even of only spin-| — there is no complete classification of entanglement. 
To be more precise, a measure of entanglement is a function on the space of states of a 
multiparticle system which is invariant under local unitary operators, i.e., unitary trans- 
formations on individual particles. Thus a complete classification of entanglement for a 
multiparticle system is a characterization of all such functions. Under the most general lo- 
cal operations assisted by classical communication (LOCC [3]), entanglement can change. 
A measure of entanglement which decreases under LOCC is called an entanglement mono- 
tone [4]. 

On two particle pure states, for example, all measures of entanglement are functions 
of the eigenvalues of the reduced density matrix (obtained by tracing the density matrix 
for the whole system over the degrees of freedom of one of the particles) , and sums of the k 
smallest eigenvalues are entanglement monotones [5] . The same information — in somewhat 
less familiar, but more algebraically convenient form — is contained in the coefficients of the 
characteristic polynomial of the reduced density matrix. These coefficients are polynomials 
in the components of the state vector and their complex conjugates. They generate the 
ring of polynomial functions invariant under the action of local unitary transformations; 
thus they completely classify two particle pure state entanglement. 

As the number of particles n increases, however, the number of independent invariants- 
measures of entanglement — grows exponentially. Complete classification rapidly becomes 
impractical. Our goal in this Letter is more modest: we seek a measure of entanglement 
which is scalable, i.e., which is defined for any number of particles; which is easily calcu- 
lated; and which provides physically relevant information. We concentrate on the case of 
spin-i particles (qubits) and begin by defining a family (parameterized by n) of functions 
on (C 2 )®". We show that each function is a measure of entanglement, vanishing exactly on 
product states. Next we evaluate this measure for several example states which illustrate 
its properties, most importantly that it measures global entanglement. This is perhaps 
best exemplified by its values on eigenstates of the antiferromagnetic Hamiltonian, for 
which we show that it is maximal only on the ground state. In a less traditional con- 
text, quantum computation relies heavily on multiparticle entangled states, particularly 
for error correction. We show that quantum error correcting code states also maximize 
our measure of entanglement. Finally, we illustrate its use in a dynamical setting, tracking 
the evolution of entanglement during a specific quantum computation. 

The Hilbert space (C 2 )® n of n qubits has a basis labelled by the 2 n n-bit strings: 
|&i . ..&„), bj E {0, 1}. For b E {0, 1}, define 

ij(b)\bi ...b n ) = Sbbjlh ...bj. ..b n ), 

where ~ denotes absence. We extend ij by linearity to be a map C 2 Cg>(C 2 )® n — > (C 2 )® n_1 . 
For u, v G (C 2 )® n-1 we can write u = ^2u x \x) and v = J2v y \y), where < x,y < 2 n ~ 1 
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are (n — l)-bit strings. Next, let 

D(U, V) = ^ \ U xVy - U y V x \ 2 , 
x<y 

the norm-squared of the wedge product of it and v. Finally, for ip G (C 2 )® 71 , define 

■ n 

Q(^) = -E £, ( t iW,»i(lW- 

n 3 = 1 

As we will see shortly, the 4/n factor provides a convenient normalization for Q. 

PROPOSITION 1. For each n G Z> 2 , Q : (C 2 )® n -> R is a measure of entanglement. 

Proof. For it, u G (C )® n_1 , D(u,v) is invariant under U(2 n_1 ). A transformation of the 
j th qubit in ^ e (C 2 )® n multiplies the j th summand in Q by the norm squared of its 
determinant. Thus each summand is invariant under TSL(2) x U(2 n_1 ), where T denotes 
the unit scalars. The intersection of these groups for all the summands is U(2) n , i.e., the 
local unitary transformations. | 

The most basic property that a measure of entanglement can have is to identify 
completely unentangled, i.e., product states. Q has this property: 

Proposition 2. Q(ip) = iff is a product state. 

Proof. Two vectors it, v G (C 2 )®"- 1 are linearly dependent iff D(u, v) = 0. Thus Q(ip) = 
implies the existence of otj G C such that ij(l)ifj = ajij(0)ip for all 1 < j < n. In particular, 

ip = |o) ®n(0)V' + \i)®n(i)ip 

= (|0) + ai|l))®n(0)i^ 
= (g®I)-(\0)®iP') 

for some g G SU(2), ip' G (C 2 )® n_1 . By Proposition 1, Q is invariant under the local 
unitary transformation g <S> I, so 

= Q(V>) = Q(|0)(g)i//) 

n 

= o + ^D(z J (o)[|o)®^],i J (i)[|o)®v1) 

J =2 

= Q(V')- 

Then, by induction, ^ is a product state. 

Conversely, if ip is a product state, then for all 1 < j < n, ij(0)ip is parallel to ij(l)ip. 
Thus Q(ip) = 0. | 



Global entanglement 



Meyer & Wallach 



Having demonstrated that Q vanishes on product states, we should now calculate it 
for some entangled states. First consider the EPR-Bohm [2,6] state (|01) — |10))/v / 2, or 
equivalently, 72 = ( 1 00) + \11})/V2. It is straightforward to calculate: 



det 1 (1 ° 



V2 V° 1 



Next, the three qubit GHZ-Mermin [7,8] state is 73 = (|000) + |lll>)/>/2- Calculating the 
invariant, we again find: 



Qhs) = 3 • I 



A* 1 f 1 

det 7f U 1 



Finally, it is now clear that for the n qubit state 7„ = (|0 . . .0) + |1 . . . l))/y/2, Q{^ n ) = 1- 
These examples demonstrate that the 4/n factor provides a natural normalization for Q. 

Proposition 3. With this normalization, < Q < 1. 

Proof. Since D(u, v) is the norm-squared of the wedge product of the two vectors u and v, 
it is bounded above by ||u|| 2 ||t>|| 2 , which takes its maximal value of \ when ||it|| 2 = ||i>|| 2 = \ 
for vectors tp = |0) <S> «i(0)V' + |1) <E> «i(l)V' with ||V'|| 2 = 1- Since there are n summands in 
Q, it is bounded above by n ■ ^ • \ = 1. | 

Proposition 2 and the calculations above show that these bounds are saturated on, 
respectively, product states and the entangled states 7 n . Of course, Q does take other 
values: Under the action of U(2) x U(2) x U(2) the Hilbert space for three qubits, C 2 ® 
C 2 (g>C 2 , decomposes into multiple orbits [9-15]. Representative states are |000) (product 
states), |0)®7 2 = |0)(|00) + |ll))/ v / 2 (and cyclic permutations), (|100) + |010) + |001))/ v / 3, 
and 73. By Proposition 2, Q(|000)) = 0. By Proposition 2 and the calculation above, 
Q( |0)(8)72) = |-| = |. A straightforward calculation gives Q(( \ 100) + |010) + |001»/ v / 3) = 
|. And we have already calculated (5(73) = 1. Thus for three qubits our measure of 
entanglement behaves in the way we would want, decreasing through states we would 
consider successively less globally entangled (and taking different values on each of these 
states, unlike the 'tangle' [16], for example, which vanishes on all but 73). In fact, on three 
qubits Q is an entanglement monotone and numerical evidence indicates that this is true 
in general [17]. 

The traditional context in which globally entangled multiparticle states occur is lat- 
tice spin systems. Consider, for example, the one dimensional spin-| Heisenberg antifer- 
romagnet [18,19] on a lattice of size n, with periodic boundary conditions, defined by the 
Hamiltonian 

n 

H n = ^ XjX j+1 + YjY j+1 + ZjZj + i, 
3=1 

where the subscripts are to be interpreted mod n, and X, Y, Z denote the Pauli matrices 
cr x , cr y , cr z , respectively. H n commutes with S z = ^2Zj, so the eigenstates of H n can be 
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labelled by their total spin S z , i.e., each eigenstate of H n is a superposition of basis vectors 
|&i . . . b n ) with \{j | bj = 1}| = s for some fixed < s < n. When s = 1, the translation 
invariance of H n implies that the eigenstates are plane waves 

1 n—l 

v4 fc) = ^5>^'|o...oio...o>, 

^ n j=0 

where the j th summand has a single 1 at the j th bit and the wave number k = 2ivm/n 
for some integer < m < n — 1. For n = 3 these plane waves are equivalent under local 
unitary transformations to the state ( 1 100) + 1 010) + (001))/^, for which we calculated 
Q = | above. In fact, for arbitrary n the entanglement of these plane waves is simply 



n ' n 2 n 2 



For s > 1 the eigenstates of H n can be computed using the Bethe Ansatz; the ground 
state has s = n/2 (for even n) [20]. Using the translation invariance of these eigenstates 
we can evaluate Q easily. The result is 

4 fn — 1\ (n — 1\ f n\ As(n — s) 



Q(S Z = s eigenstate of H n ) = ^ • n ■ ^ ^ ^ ^ 



n 2 



Notice that for the ground state this entanglement measure is maximal, Q = 1. This 
result contrasts with O'Connor and Wootters' calculations of the 'concurrence' C in these 
states [21]: they find that C is not maximal on the ground state, but rather for s/n ~ 0.3 
as n — > oo. This difference is due to the fact that C is really a measure of two particle 
entanglement, even when generalized to multiparticle states, while Q is a global measure 
of multiparticle entanglement. 

Highly entangled multiparticle states also occur in the relatively new context of quan- 
tum error correcting codes. In fact, the code subspace for an additive code can be described 
as the space of ground states of the Hamiltonian formed by the sum of the stabilizers [22] . 
For example, the code subspace for a 5 qubit code [5,1,3] encoding 1 qubit against single 
bit errors [23,24] is the space of ground states of the translation invariant Hamiltonian on 
a one dimensional lattice of five qubits: 

5 

#[5,1,3] = XjZ j+1 Zj +2 x j+?n 

3 = 1 

where the subscripts are to be interpreted mod 5. The space of ground states is two 
dimensional — which is why it can encode 1 qubit. A basis is 

|0) i-> [|00000) - (|11000) + eye.) + (|10100) + eye.) - (|11110) + cyc.)]/4 
|1) i-> [|11111> - flOOlll) + eye.) + (|01011) + eye.) - (|00001) + cyc.)]/4, 



5 



Global entanglement 



Meyer & Wallach 



where "eye." indicates cyclic permutations. From these equations it is straightforward 
to calculate that Q = 1 for all states in this code space. Here the difference from the 
concurrence is even more dramatic: C vanishes on the code subspace since tracing over all 
but two qubits leaves a reduced density matrix proportional to the identity [25]. 

Shor's original 9 qubit code protecting 1 qubit against single qubit errors [26] can 
also be described as the ground state subspace of a lattice Hamiltonian — for a lattice 
triangulating M.P 2 [27]. In this case a basis for the code space is 

|0> i-> (|000) + |lll))® 3 /3 
|1> i-> (|000> - |lll»® 3 /3. 

Calculating Q for the states in this subspace we find that again it is maximal, despite 
the fact that these states decompose into products of three qubit factors. So Q does not 
distinguish all sub-global entanglements; this is a consequence of using a single invariant. 
Finer resolution requires a more complete set of invariants, and in general, higher degree 
polynomials [9,10,12,15,28]. 



Nevertheless, as we have seen, Q pro- 
vides useful information about global en- 
tanglement in certain contexts. Further- 
more, in dynamical problems, Q quanti- 
fies the evolution of entanglement. Con- 
sider Grover's algorithm [29] , for example: 
Given a,j e {0, 1}, 1 < j < n, define 
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Figure 1. Entanglement in Grover's algorithm 
for 10 qubits as a function of number of iterations. 



and then extend U a by linearity to a map 
(C 2 )® n -> (C 2 )® n . The goal of Grover's 
algorithm is to convert an initial state of 

n qubits, say |0. . .0), to a state with probability bounded above \ of being in the state 
\a\ . . .a n ), using U a the fewest times possible. Grover showed that it can be done with 
0(\/2™) uses of U a by preparing the state 

2™— 1 

-L V \x) =H® n \0...0) where H=^=(\ \ V 

and then iterating the transformation H® n U H® n U a on this state [29]. The initial state is 
a product state, as is the target state \a\ . . . a n ), but intermediate states ip(k) are entangled 
for k > iterations. We can evaluate Q on these states to quantify this entanglement: 



nfinW a( N A cos *Vkf. a costf fc y 
Q(W)=4( y -lj 17 - T (s 1 n^-- 7 _j , 



cos 2 9k 



cos 61 
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where 9k = (2k + 1) esc 1 (\/N) and N = 2 n . The results are plotted in Figure 1 for 
n = 10: the entanglement oscillates, first returning to close to at 

jfc= [If 

^2 V2csc- 1 (v / iV) 

where [•] denotes 'closest integer to'; this is when the probability of measuring \a\ . . .a n ) 
is first close to 1. 

Three qubit states, eigenstates of lattice Hamiltonians, quantum error correcting code 
subspaces, and the intermediate states in Grover's algorithm all illustrate how a measure 
of multiparticle entanglement such as Q provides insight into global properties of quantum 
multiparticle systems. While Q has the satisfactory properties of Propositions 2 and 3, is an 
entanglement monotone on three qubits, and is a straightforwardly computable polynomial, 
it is in no sense a unique measure of multiparticle entanglement. A more complete (but 
still partial) characterization can be obtained by also using some of the other measures 
which have been proposed, like the concurrence [21], the closely related n-tangle [30], 
the Schmidt rank [31], the negativity [32], etc. Each emphasizes a specific feature of 
multiparticle entanglement and describes a different physical property. We anticipate that 
multiparticle entanglement measures — whose current development is largely motivated by 
quantum computation — will contribute to the understanding of the physics of quantum 
multiparticle systems more generally [33-35]. 
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